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Abstract 

We demonstrate that, at a 500 GeV e + e~ collider, inclusion of the ZZ- 
fusion process for production of a light standard-model-like Higgs boson can 
substantially increase the precision with which the ZZh coupling can be de- 
C^ - ■ termined (using the model-independent recoil mass technique) as compared 

^ \ to employing only Zh associated production. 

o 

2£ ' 1 Introduction 

Once a neutral Higgs boson (h) is discovered, determining its coupling to Z bosons 
(9zzh) i s OI fundamental importance. [] It is this coupling which most directly re- 
flects the role of the h in electroweak symmetry breaking. In the minimal Standard 
Model (SM), where the Higgs sector consists of a single Higgs doublet field, there 
is only one physical Higgs boson eigenstate, with coupling g ZZh = gmz/ cosOw, 
where g is the SU(2) coupling and 9w is the weak mixing angle. In contrast, in 
a theory with many scalar doublets and/or singlets, the ZZ couplings of the in- 
dividual neutral Higgs bosons (hi) are generally reduced in magnitude, but must 
obey the sum rule Y^i9zzhi = [9 m z/ cos^] 2 . The sum rule becomes still more 
complicated if triplet Higgs representations are present. Precise determination of 
9zzh f° r each and every observed h will therefore be crucial to knowing whether or 
not we have found all the Higgs bosons that participate in electroweak symmetry 
breaking and to understanding the full structure of the Higgs sector. 

In e + e~ collisions, the dominant Higgs boson production diagrams involving 
the ZZh coupling are of two types: 

e + e~ -> Zh (1) 
e + e~ -> e + e~Z*Z* -> e + e"/i. (2) 



1 For recent reviews of Higgs boson phenomenology, see e.g., Refs. |^j. 
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There is (constructive) interference of the amplitude for e + e~ — > Zh — ► e + e~h with 
that for (Q). However, it is desirable both for simplicity and in order to maximize 
experimental accuracy for the g ZZh determination to impose cuts such that this 
interference is very small; the Zh (JED and ZZ-fusion (H) amplitudes can then be 
considered as leading to effectively independent production processes. f\ The main 
point of this paper is to demonstrate that, in the case of a light Higgs boson, when 
the e + e~ collider is operated at full energy (e.g. \/s = 500 GeV), the ZZ-fusion 
production mode will make a crucial contribution to the determination of the ZZh 
coupling. Indeed, by combining the ZZ-fusion and Zh production processes we 
find that the error for g ZZh achieved at \/s = 500 GeV can be competitive with 
that which is attained using Zh associated production (alone) at lower y/s near 
the maximum in the Zh cross section. Further, if the linear collider is run in the 
e^e - mode, the only source of h production, and only means for determining the 
ZZh coupling is from ZZ-fusion e - e~ — ► e~e~h |J. 

Processes (0) and (fj) have quite different characteristics. For a lighter Higgs 
boson and lower center of mass energy (y/s) process P) dominates (with a maximal 
cross section at y/s ~ m z + v2 nth) , while for a heavier Higgs boson or higher y/s, 
process (|2|) becomes more important (the cross section increasing logarithmically 
with \/s). This is shown in Figure [I], where we present the total cross section 
for e + e~ — > e + e~h (solid) as a function of y/s for several Higgs boson masses, 
rrih = 80, 120, 160 and 200 GeV. Dashed curves present the contribution only from 
e + e~ — > Zh with Z — > e + e~. Note that the ZZ fusion cross section becomes larger 
than that of Zh for y/s > 300 GeV. 

In both the Zh and ZZ-fusion channels, the Higgs signal can be easily detected 
for nth ^ (0.7 — 0.8) y/s by reconstructing the Higgs mass peak via the h decay 
products. However, in order to determine the ZZh coupling in a model-independent 
manner (i.e. independent of the h's branching ratio to any particular channel), it 
will be crucial to identify the Higgs signal through the "recoil mass" variable, 

M r 2 ec = s + Mf e — 2\/s (E t+ +E e -), (3) 

where M« is the invariant mass of the final state lepton pair and the Ei are the 
lepton energies in the cm. frame; here, I = e, [i is possible for process (0) [] while 
only t = e is relevant for process (0). Due to detector resolution effects, M rec will 

2 Nonetheless, our calculations will always employ the full SM matrix elements j|, including all 
interfering diagrams, for any particular signal final state. Full matrix elements are also employed 
for background processes contributing to any particular final state. However, interference between 
the signal and background is neglected; this is an excellent approximation when considering a 
very narrow light Higgs boson. 

3 One could also consider reconstructing an M rec peak using the Z — > qq (hadronic) decays 
in Eq. (^) in order to increase the signal statistics. However, the energy /momentum resolution 
for jets is much worse than for leptons and the backgrounds in the hadronic decay channels arc 
larger, implying a less sharp signal peak above background. Thus, we will consider only the 
leptonic modes. 
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Figure 1: Total cross section for e + e~ — > e + e~h (solid) as a function of y/s for 
rrih = 80, 120, 160 and 200 GeV. Dashed curves present the contribution only from 
e + e~ — ► Zh with Z — > e + e~. 

display a peaked distribution of finite width, much broader than the physical width 
of a light Higgs boson, centered about m^. If we can measure the inclusive cross 
section associated with the M rec peak in such a way that there is small sensitivity 
to the h decay, then we can obtain a direct determination for the ZZh coupling 

9zzh- 

Since we allow the Higgs to decay to anything, our background is composed of 
many processes. For i = e, for example, we must consider all processes of the type: 

e + e- -> e + e"X, (4) 

with X = £ + £~ , r + r _ , vv and qq. 

Many analyses of the process ([!]) in the inclusive M rec context have appeared in 
the literature; see, for example, || 0, 0, || and references therein. However, pro- 
cess (H) has received limited attention p], [10| . In particular, complete background 



computations for the inclusive signal are given for the first time in the present 
paper. 



3 



The rest of the paper is organized as follows. In Sec. 2 we explore the relative 
importance of the processes (JT|) and (||) for determing g 2 ZZh (assuming a SM-like 
h), and comment on the e~e~ — > e~e~h ZZ-iusion production mode. Section 3 
summarizes our results and their implications for the relative importance (for the 
9zzh determination) of running at y/s = 500 GeV vs. lowering ^/s to a value near 
the maximum in the Zh cross section. 



2 Determining the Coupling g 2 zzh 

In order to reconstruct the signal peak in M rec via Eq. (|3]), we must assure that the 
charged leptons are detected. Thus, we impose the following "basic" acceptance 
cuts 

| cos 0*| < 0.989, VriP) > 15 GeV, p T (t + t~) > 30 GeV, 

M u > 40 GeV, and M rec > 70 GeV. (5) 

The polar angle cut roughly simulates the detector acceptance for a beam hole 
of 150 mrad ||, and the cut on Mu is imposed in order to suppress events from 
photon conversion. 

The sharpness of the reconstructed M rec peak is determined by the momen- 
tum/energy resolution for the charged leptons. The energy of an electron (but 
not that of a muon) can be measured in the electromagnetic calorimeter. The 
momentum of either a muon or an electron can be determined from a tracking 
measurement of its curvature in the magnetic field of the detector. We consider 
two possibilities for the energy resolution of the electromagnetic calorimeter: 



I: NLC/EM AE/E = 12%/y/E © 1%; 

II: CMS/EM AE/E = 2%/\[E~ © 0.5%. 

where © denotes the sum in quadrature and E is in GeV. Case I is that currently 
discussed for the NLC electromagnetic calorimeter B|; case II is that of the CMS 



lead-tungstate crystal [11|. We also consider two possibilities for the momentum 



resolution from tracking: 

III: NLC/tracking Ap/p = 2 x 10~ 4 p © 1.5 x lO^V^sm 2 0; 

IV: SJLC/tracking Ap/p = 5 x 10 _5 p © 10~ 3 , 
with p in GeV. Case III is that specified for the typical NLC detector in [H] and 



case IV is that quoted for the "super- JLC" (SJLC) detector design |T2 . 

It is important to note that for an electron the tracking determination of its 
momentum is not statistically independent of the electromagnetic calorimeter mea- 
surement of its energy. Thus, the two measurements cannot be combined; one 
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should use whichever measurement provides the best result. In order to provide a 
clean comparison of the different possibilities, we will present results in which we 
analyze all events using either the calorimeter energy measurement or the tracking 
measurement; i.e. we do not choose the best measurement on an event-by-event 
basis. 

It is useful to compare the fractional resolutions, r = AE/E (E ~ p for I = 
e,/i), for the different cases to one another as a function of energy/momentum. 
Using 9 = 90° in III, one finds that rni > T\ for any E > 70 GeV (lower E 
for 9 in the forward or backward direction) and that rj > rrv > rn for E > 
100 GeV. At yfs = 500 GeV, lepton energies above 100 GeV are typical for the light 
Higgs masses studied here; then, NLC/tracking will not be useful for an electron, 
whereas SJLC/tracking might be, depending upon the EM calorimeter. Since 
the muon energy /momentum can only be measured by tracking, NLC/tracking 
will result in the Zh — > fi + fi~h channel having a poorer signal to background 
ratio than either the Zh — ► e + e~h or the ZZ-fusion e + e~h channel. On the 
other hand, if the machine energy is lower, e.g. near the peak in the Zh cross 
section for small m^, electron energies are smaller, and for the majority of events 
the NLC/tracking measurement of the electron energy is competitive with the 
NLC/EM measurement. This will be apparent from the figures in the next section. 

To suppress the SM background more effectively and for physics clarity, it is 
beneficial to divide the study into the two natural categories: the Zh associated 
production and ZZ fusion processes. 

2.1 e + e~ — ► Zh — ► e + e~h and ^L + fj,~h 

We first discuss a linear collider with y/s = 500 GeV. In order to isolate the Zh 
class of events we require 

\M u -m z \ < 10 GeV , |cos0*|<O.8. (6) 

The mass cut largely eliminates the leading background from e + e~ — > W + W~ and 
the polar angle cut helps reduce the other large background from e + e _ — > ZZ. 
Since the Z boson in the signal is not only central, but also energetic with E z = 
(s — ml + m 2 z )/2y/s ~ y/s/2, we can further reduce the background by imposing 
a cut of 

p T (e + e-) > 80 GeV . (7) 

The signal to background ratio is improved after the cuts to the extent that the 
Higgs must be nearly degenerate with the Z for the M rec peak to have any signif- 
icant background. Figure |] presents the recoil mass distributions, da/dM rec , for 
e + e~ — > e + e~X at y/s = 500 GeV. Results for different energy /momentum reso- 
lutions, cases I-IV, are shown in the four different panels. The solid and dashed 
curves give the Zh signal for = 90 and 120 GeV, respectively. The dotted line 
is the summed SM background. 
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Figure 2: Recoil mass distributions for e + e~ — > e + e~X at \fs = 500 GeV. The 
solid and dashed curves give the Zh signal for m h = 90 and 120 GeV, respec- 
tively. The dotted line is the summed SM background. Results for the different 
energy /momentum resolution, cases I- IV, are shown in the four different panels. 
The cuts of Eqs. (J5|), (@) and (0) have been imposed. 

The e + e _ — > Zh cross section reaches a maximum near y/~s ~ mz + V% n^h- A 
relevant question is how much improvement is possible by running the machine at 
a lower energy nearer the maximum cross section, as would be possible once is 
known (either from LHC or NLC data). To illustrate, we consider y^i = 250 GeV. 
Figure [| shows the recoil mass distributions, da/dM rec , similar to those in Fig. |2] 
but for yfs = 250 GeV. We see not only that the cross section rate is larger, but 
also that the signal peak is much sharper because of the better determination for 
lepton energy/momentum at this lower y/s. 

We estimate the relative statistical error for the cross section measurement as 

R = VSTB/S , (8) 

where 5* and B are the numbers of signal and background events for a given lu- 
minosity; neglecting systematic uncertainty from correcting for our cuts, this is 
also the error on the coupling g 2 ZZh . In Table [l] we compare the errors (as a func- 
tion of rrih) in the Zh mode, with Z — > e + e~, for different resolution choices, 
taking ^/s = 500 GeV and 250 GeV and assuming an integrated luminosity of 
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Figure 3: Recoil mass distributions for e + e~ — > e + e~X at y/s = 250 GeV. The 
solid and dashed curves give the Zh signal for m h = 90 and 120 GeV, respec- 
tively. The dotted line is the summed SM background. Results for the different 
energy /momentum resolution, cases I- IV, are shown in the four different panels. 
The cuts of Eqs. (|5|) and (|J) have been imposed. 

L = 200 fb -1 . At = 500 GeV, the error ranges from ~ 15% at m h = 80 GeV 
to ~ 7% at m h = 140 GeV. In general, the accuracy is not greatly affected by the 
detector. The exception is case III which yields poorer results than the other cases. 
Results for the g 2 ZZh error at yfs = 250 GeV are ~ 6% (more or less independent of 
resolution choice) for most values of m^, worsening to ~ 7% at the Z peak. Thus, 
by running the machine at an energy near the maximum in the Zh cross section 
one should be able to improve the accuracy for this particular mode by about a 
factor of 2. 

We compute the net error for the g\zh determination in the Zh channel by 
including measurements for both the Z —>■ e + e~ and the Z — > final states. 

The net error (_R ne t) is given by 

R net = R e 2 + ^ ■ ( 9 ) 

In obtaining the error for the e + e~ final state, we compare results found using the 
EM calorimeter to those found using tracking and adopt the superior choice. Thus, 
for example, if we assume NLC/EM and NLC/tracking, this means using case I 
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Table 1: Percentage accuracy for g 2 ZZh based on Zh channel e + e~ — > e + e~h cross 
section measurements with L = 200 fb _1 assuming (a) y/s = 500 GeV and (b) 
y/s = 250 GeV. Results for the four different lepton energy/momentum resolutions 
are shown. The cuts of Eqs. (^|) and (||) have been imposed for both y/s = 500 GeV 
and y/s = 250 GeV. For y/s = 500 GeV we have imposed the additional cut (|7]). 
Note that only one kind of lepton is counted here. 
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resolution 


mass bin 


m h (GeV) 
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I 


(m h ± 10) 


15% 


15% 
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8.3% 


7.2% 


II 


(m h ± 10) 


11% 


12% 


10% 


6.2% 


6.4% 


III 


(m h ± 10) 


19% 
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15% 
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IV 


{m h ± 10) 


11% 


12% 
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6.8% 


6.5% 




(b) 


yfs = 250 GeV 






I 


(m h ± 10) 


5.7% 


6.7% 


6.3% 


4.8% 


6.6% 


II 


(m h ± 10) 


5.8% 


6.8% 


6.7% 


4.8% 


6.5% 


III 


K ± 10) 


5.6% 


6.7% 


6.7% 


4.7% 


6.7% 


IV 


(m h ± 10) 


5.8% 


6.7% 


6.7% 


4.7% 


6.7% 



Table 2: The percentage accuracy for g ZZh obtained by combining [via Eq. (|])] 
results for the Zh cross section measurement in the e + e~ — > e + e~h (NLC/EM) 
channel and the e + e~ — > fi + n~h (NLC/tracking) channel, assuming L = 200 fb _1 
and y/s = 500 GeV or y/s = 250 GeV. The cuts of of Eqs. (Q) and (|6]) have been 
imposed for both energies and for y/s = 500 GeV we impose the additional cut (]?]). 
The mass bins of Table have been employed. 



GeV) 
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10% 
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4.7% 


4.6% 


3.4% 


4.7% 



s 



results for the Z — > e + e~ final state and case III results for the Z — > final 
state. The results are illustrated in Table At y/s = 250 GeV our results indicate 
that an error for g% Zh of less than 5% can be expected. [] 

2.2 e + e~ — > e + e~h and e~e~ — ► e~e~h via Z Z-fusion 

The major advantage for the ZZ-fusion channel (|2]) over the Zh channel ([!]) is that 
the cross section increases logarithmically with energy; at yfs = 500 GeV and for 
rrih = 120 GeV it is about 10 fb as compared to the Zh cross section of about 2.5 
fb. To remove the large ZZ background, we require, in addition to the basic cuts 
ofEq. (D, 

M ee >100GeV. (10) 

The only penalty from the M ee cut is a 20% decrease in the signal rate due to 
elimination of the constructive interference of the ZZ-fusion amplitude with the 
Zh amplitude. 

Figure £| presents signal and background curves after the cuts of Eqs. @ and 
(|T0"D. Although some background persists, f\ the much larger signal rate makes 
a good measurement of the cross section feasible. The corresponding results for 
the error of the g ZZ f l determination are presented in Table [3|. Using NLC/EM 
calorimetry, case I, the error ranges from 10% to 7%. As in the Zh channel, there 
is little difference between results for the resolution cases I, II and IV, whereas case 
III yields substantially poorer results. 

The ultimate accuracy for g 2 ZZh at y/s = 500 GeV is obtained by combining the 
Zh and ZZ-fusion channel results. In Table f|, we present the error as a function 
of rrih, assuming that NLC/EM energy resolution is employed for the e + e~h final 
states of the Zh and ZZ-fusion channels and that NLC/tracking resolution is 
employed for the n + n~h final state of the Zh channel. Note that the achievable 
accuracy is competitive with that obtained at y/s = 250 GeV using the Zh mode 
alone (the ZZ-fusion mode being not useful at this low energy), especially at the 
larger values. 

At an e~e~ collider, Higgs production is entirely from the ZZ-fusion process, 
er er — ► e~e~h. The results presented in Table || are essentially applicable, except 
that the background level is slightly smaller in the e~e~ case [H. 



4 This, and other results obtained here for the Zh mode are generally consistent with those 
of Refs. (|, 0, ||, |9| for mh <; 110 GeV, when the same electromagnetic-calorimeter/tracking 
resolution assumptions are made. For lower mh, we find higher backgrounds as compared to the 
estimates made in Ref. resulting in larger errors. 

5 The signal-to-background ratio can be further enhanced by stronger cuts, but the signal rate 
is also reduced and the error in the determination of g 2 zzh is not improved. 
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Table 3: Percentage accuracy for g 2 ZZh based on the ZZ-fusion channel e + e — > 
e + e~h cross section measurements with L = 200 ftT 1 at y/s = 500 GeV. Results for 
NLC/EM, CMS/EM and SJLC/tracking (resolution cases I, II and IV) are shown; 
NLC/tracking yields much poorer results than NLC/EM. The cuts of Eqs. (|), and 
fllUD have been imposed. 
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Table 4: Combined percentage accuracy for g ZZh as obtained by including the 
three measurements: e + e~ — > e + e~h for the Zh and ZZ-fusion (NLC/EM) plus 
e+e~ -> fM+ii-h for the Zh (NLC/tracking) at ^ = 500 GeV and L = 200 ftT 1 . 
Cuts as specified in Tables p] and |] are imposed. 
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Figure 4: Recoil mass distributions for e + e~ — > e + e~X at = 500 GeV. Results 
for m h = 90 GeV (solid), and 120 GeV (dashed) are shown. Results for the different 
energy /momentum resolution, cases I- IV, are shown in the four different panels. 
We impose the cuts of Eqs. (|S|) and (fTTH). 

3 Summary and Conclusions 

We have investigated the precision with which the Higgs boson to ZZ coupling 
(9zzh) can be determined by employing the recoil mass distribution in £ + £~h (i = 
e, /i) final states coming from the Zh and ZZ-fusion production channels at an 
e + e _ collider or ZZ-fusion (alone) at an e~e~ collider. From the (fully inclusive) 
recoil mass distribution a direct determination of g 2 ZZh , that is independent of 
any assumptions regarding the branching ratio for the Higgs boson to decay into 
any particular channel, can be made. We considered a number of electromagnetic 
calorimetry and tracking resolution options for the detector. Results for the four 
options considered are similar, except that the tracking specified for the 'typical' 
NLC detector yields poor results at y/s = 500 GeV. For the calorimetry and 
tracking resolutions specified for this typical NLC detector, we find the following 
results for the error of the g\zh determination for mh in the range 80 — 140 GeV: 

• If the e + e~ collider is run at a/s = 500 GeV, the ZZ-fusion production mode 
yields smaller statistical error than does the Zh production mode (even after 
combining both e + e~h and fi + fi~h final states in the latter case). 
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• Taking both Zh channel and ZZ-fusion channel measurements at y/s = 
500 GeV and 200 fb -1 into account, the combined accuracy ranges from 7.5% 
at m h = 80 GeV to 4.7% at m h = 140 GeV (Table |). 

• Running at a lower energy near the Zh cross section maximum is possibly 
fruitful. For accumulated luminosity of 200 ftT 1 at y/s = 250 GeV, the 
error found using the (dominant) Zh production mode is less than 5% for 
80 < m h < 140 GeV (Table §. 

• The statistical error of a measurement in the ZZ-fusion e~e~ — > e~e~h 
channel would be somewhat better than that for a measurement in the ZZ- 
fusion e + e~ — > e + e~h channel (Table ^), due to smaller backgrounds. 

The most significant implication of our results is that it may not be necessary, 
or even appropriate, to run at low energy in order to obtain the best possible 
accuracy for the g 2 ZZh determination. When running at the full energy of y/s = 
500 GeV, if m h > 120 GeV then (for NLC/EM and NLC/tracking) the combined 
Zh and ZZ-fusion error of Table ^ is very close to that obtained when running at 
t/s = 250 GeV for the same integrated luminosity. Further, it seems likely that it 
will prove desirable from the point of view of other physics to accumulate more 



luminosity at ^fs = 500 GeV than at ^ = 250 GeV, and it is also possible [13 
that the instantaneous luminosity at the lower energy will be lower (assuming that 
the interaction region is optimized initially for -^i = 500 GeV). In either case, the 
lower energy running might not significantly improve the g\zh determination even 
if m h < 120 GeV. 

Thus, we conclude that use of the ZZ-fusion mode (as well as the Zh mode) 
should provide a very valuable increase in the accuracy that can be achieved for 
the determination of the ZZ coupling of a SM-like Higgs boson at a lepton collider 
operating at high energy. 
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